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Abstract — The spin torque assisted thermal switching of the 
single free layer was studied theoretically. Based on the rate 
equation, we derived the theoretical formulas of the most likely 
and mean switching currents of the sweep current assisted 
magnetization switching, and found that the value of the exponent 
b in the switching rate formula significantly affects the estimation 
of the retention time of magnetic random access memory. Based 
on the Fokker-Planck approach, we also showed that the value 
of b should be two, not unity as argued in the previous works. 

Index Terms — spintronics, thermal stability, Fokker-Planck 
equation, theory 



I. Introduction 

MAGNETIC random access memory (MRAM) using 
tunneling magnetoresistance (TMR) effect [1],[2] and 
spin torque switching [3], [4] has attracted much attention 
for spintronics device applications due to its non-volatility 
and fast writing time with a low switching current. A high 
thermal stability (Ao) (more than 60) of magnetic tunnel 
junctions (MTJs) is also important to keep the information 
in MRAM more than ten years. Recently Hayakawa et al. 
[5] and Yakata et al. [6], [7] respectively reported that the 
anti-ferromagnetically (AF) and ferromagnetically (F) coupled 
synthetic free layers show high thermal stabilities (Ao > 80 
for AF coupled layer and Ao = 146 for F coupled layer) 
compared to a single free layer. 

The thermal stability has been determined by measuring the 
spin torque assisted thermal switching of the free layer and 
analyzing the time evolution of the switching probability by 
Brown's formula [8] with the spin torque term. The theoretical 
formula of the switching probability is generally given by 
P = 1 - exp[- $1 dt'v(t% where u(t) = / exp[-A (l - 
I/I c ) b ). Here, /o, /, and I c are the attempt frequency, current 
magnitude, and critical current of the spin torque switching at 
zero temperature, respectively, b is the exponent of the current 
term in the switching rate v, and was argued to be unity by 
Koch et al. in 2004 [9]. On the other hand, recently, Suzuki et 
al. [10] and we [11], [12] independently studied the spin torque 
assisted thermal switching theoretically, and showed that the 
exponent b should be two. Since the estimation of the thermal 
stability strongly depends on the value of 6, as discussed in this 
paper, the determination of b is important for the spintronics 
applications. 

In this paper, we study the spin torque assisted thermal 
switching of the single free layer theoretically. In Sec. HU we 
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derive the theoretical formulas of the most likely and mean 
switching currents of the sweep current assisted magnetization 
switching, and study the effect of the value of the exponent 
b on the estimation of the retention time of the MRAM. In 
Sec. HID the differences of the theories in Refs. [9], [10], [11] 
are discussed by analyzing the solution of the Fokker-Planck 
equation. Section [IV] is devoted to the conclusions. 



II. Theory of Magnetization Switching due to 
Sweep Current 

In this section, we consider the spin torque assisted thermal 
switching of the uniaxially anisotropic free layer, which has 
two minima of its magnetic energy. At the initial time t = 
0, the system stays one minimum. From t — 0, the electric 
current I(t) = xt is applied to the free layer which exerts the 
spin torque on the magnetization and assists its switching. In 
this section, the current is assumed to increase linearly in time 
with the sweep rate xr, as done in the experiments [7], [13], [14]. 
The magnitude of the current I(t) = xt should be less than 
I c because we are interested in the thermally activated region. 
The time evolution of the survival probability of the initial 
state, R(t), is described by the rate equation, 



dR(t) 
dt 



-v(t)R(t), 



(1) 



where the switching rate v(t) is given by 



v(t) = /oexp 



-An 1 - 



I(t) 



(2) 



We assume that the attempt frequency is constant, b is the 
exponent of the current term, (1 — I/I c ). The switching 
probability is given by P(t) = 1 — R(t). Also, we define 
the probability density p(t) by p(t) = -dR/dt = dP/dt. 
Equation ([]]) describes the escape from one equilibrium to the 
others in many physical systems, and the value of b reflects 
their energy landscape: 6=1 for the Bell's approximation 
[15], b — 3/2 for the linear-cubic potential [16], and 6 = 2 
for the parabolic potential [17], [18]. The determination of the 
value of b has been discussed not only in spintronics but also 
the other fields of physics [19]. The form of Eq. (ffj is the 
special case of the model of Garg (a in Ref. [20] corresponds 
to 1 - 6). 



2 



(a) 

£1.0 



| 0.6 
LO-4 

I) 

| 0.2 

i 

; 0, 



b=2 b=1 



0.2 0.4 0.6 0.8 1.0 
time (s) 




0.2 0.4 0.6 0.8 
time (s) 



Fig. 1. The time evolutions of (a) the switching probability P(t) and (b) 
its density p(t) for 6 = 1 (solid) and 6 = 2 (dotted). 

The solution of Eq. (Q} with the initial condition R(0) = 1 
is given by 



R(t) =exp • 



bxA 



1/6 




(3) 

where 7 (/J, 2) = Jl dtt^~ 1 e~* is the lower incomplete T 
function. Figure Q] the time evolutions of (a) the switching 
probability P(t) and (b) its density p(t). The values of the 
parameters are taken to be fa = 1.0 GHz, I c = 1.0 mA, 
x = 1.0 mA/s, and Ao = 60, respectively, which are typical 
values found in the experiments [6], [7], [13], [14]. As shown, 
P(t) suddenly changes from to 1 at a certain time t = t 
at which p(t) takes its maximum. We call t the switching 
time. The switching time i is determined by the condition 
(dp(t)/dt) t -i — 0, i.e., dv/dt — v 2 , and is given by 



xt 1 
P= 1 -A- l ° S 



xA 



(4) 



for 6 = 1, and 



T = 



6-1 
"6A7 



plog 



1 



bxA 



b/(b-i) 



1/6 



i/b 



(5) 

for 6 > 1. Here plog(z) is the product logarithm which 
satisfies plog(z) exp[plog(z)] = z. For a large z ^> 1, 
plog(z) ~ log 2:, and t (6 > 1) can be approximated to 



>rf 1 , 



1 



1-1/6 



6^A 



1/6 
j 



(6) 



The current at t = t, I(t) = xt, is the most likely switching 
current for the thermal switching. Since we are interested in 
the switching after the injection of the current at t = 0, t 
should be larger than zero. Thus, the above formula is valid 
in the sweep rate range x > x c , where the critical sweep rate 
x c is given by 

fnlc -An 



6An 



(7) 



The value of x c estimated by using the above parameter values 
is on the order of 10~ 19 mA/s, which is much smaller than the 
experimental values (0.01 — 1.0 mA/s in Ref. [14]). Thus, the 
above analysis is applicable to the conventional experiments. 
We also define the mean switching current (J) by 
,.1 



{!) = 



dRI = 



' , dR 
dt——xt = x 
dt 



dtR. (8) 



Since p(t) takes its maximum at t = t, we approximate that 
dv 

t=t 



v(t) + 



dt 



(t - t) = v [1 + v (t - t)] ^ ve*^, 

(9) 



where v — Then, R(t) = exp[— J Q dt'v{Jt')\ can be 

approximated to 



where A 



R(t) ~ exp{-A [exp (vt) - 1]} 
~ vt . Thus, (/) is given by 



(/) ~ xc A j dtexp (-Ae m ) = ^-^(A), 



(10) 



(11) 



where Ep(z) — f^dte zt /t 13 is the exponential integral. It 
should be noted that E\(A) is expanded as [21] 



B 1 (A) = -7-logA-^ 



fe=i 



(~A) fc 
kk\ ' 



(12) 



where 7 = 0.57721... is the Euler constant. In general, the 
moment (7 n ) = dRI n = nx n J °° dtRt 71 ^ 1 is given by 



nx n e A 







■ x\ 
n ( — I e' 

. v . 



dt t n ~ l exp (~Ae l ' t ) 



(13) 



Then, the standard deviation of the current, 07 = y/ (I 2 ) — (I} 2 , 
is given by 
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2+7(2- 7 ) + logA(2-2 7 -logA) 
+7(7-37) + logA(7-67-31ogA) 



^(A 4 )} 



(14) 



Since the thermal stability can be estimated by evaluating 
the parameter v, as shown below, let us derive the relations 
between v and experimentally measurable variables. The dif- 
ference between the most likely switching current = xi 
and mean switching current (I) is given by 



</>-/(*) = -^T 



7 + E 



(-A)' 



fc=i 



fcfc! 



For 6 = 1, v and i>i are, respectively, given by 

>fA 



ljxi. (15) 



(16) 



Vt = An 



, 1 1 



A V ^ A o 



(17) 



As shown in Refs. [22], [23] I(t)/I c — vt//S. is around 
0.4 ~ 1.0 in the experimentally reasonable temperature and 
sweep rate regions (so called fast pulling regime or Garg's 
limit [24], [25]). Thus, we can approximate that A = eT vt ~ 



3 



and e = e c ~ e c 
for b = 1 is given by 



(I)-I(t) = -7 



1 for A > 1. Then, (I) - I(t) 

Ic. 



A " 



(18) 



Similarly, for b > 1, by using the approximation plog(z) 
log z, v and vt are, respectively, given by 



6-1 



-, \ l-l/b 



1 - 1/b bxA 1/b 



o 



(19) 
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-1 


i log 







6-1 



1-1/6 



Zoic 



6><A 



1/6 



1/6N 



6-1 



1-1/6 



1/6 



Then, (7) — I(t) for 6 > 1 is given by 



6 V" 1 ' 6 J B 



6-1 



6A 



1/6- 



(20) 



(21) 



[(/) — J(i)]/7 C is approximately zero for a sufficiently high 
thermal stability (Ao 3> 1) which means a narrow width of 
the probability density. We also find 



v [{I)-I(t)] 



-7 = -0.57721... 



= 1.28254... 



(22) 



(23) 



for arbitrary 6 and Ao 3> 1. We numerically verify Eqs. 
and ( 1231 among the temperature region < T < 500 K, where 
the values of the parameters are same with those in Fig. [TJ 
(A oc l/T is taken to be 60 for T = 300 K). Equation (f22| 
or d23l can be used to determine the value of v experimentally. 
Otherwise, v can be estimated by using the relation 



1 dR 

R~dJ 



dt 



logi? 



(24) 



Let us discuss the effect of the value of 6 on the estimation 
of the retention time of MRAM. We assume that the value of I c 
is experimentally determined by some other experiments [5], 
Then, the unknown parameter in Eq. ([TBI or ([19) is only the 
thermal stability. As mentioned above, v can be experimentally 
determined by using Eq. (l22l . (1231 . or (l24l . By setting v(b = 
1) = v{b = 2), we found that the estimated values of the 
thermal stability with 6 = 1 (Ai) and 6 = 2 (A2) satisfy 
the relation Ai = ^/TK^. Let us define the retention time 
of MRAM by t* = e A °//o- Then, the ratio of the estimated 
values of the retention time by b = 1 (tj) and 6 = 2 (t^) is 
given by fjjj/i* = e A2 ~v / 2£i" j which is on the order of 10 21 
for A 2 = 60 and increases with increasing A 2 . Thus, the 
determination of the value of 6 is important for the accurate 
estimation of the retention time of MRAM. 



III. Comparison with Theory of Koch et al. 

In this section, we investigate the difference of the value of 6 
between Koch et al. [9] and Refs. [10], [11], [12] by comparing 
the solutions of the Fokker-Planck equation, and show that 
6 should be two. For simplicity, in this section, the current 
magnitude is assumed to be constant in time [9],[10],[1 1],[12]. 

First of all, it should be mentioned that the analytical 
solution of the switching probability can be obtained only for 
the two special cases. The first one is the uniaxially anisotropic 
system [10]. The second one is the in-plane magnetized thin 
film in which the switching path in the thermally activated 
region is completely limited to the film plane, and thus, the 
effect of the demagnetization field normal to film plane is 
neglected [11]. In these systems, the magnetization dynamics 
can be described by one variable (the angle from the easy 
axis, 9), although, in general, the magnetization dynamics is 
described by two angles (the zenith angle 9 and azimuth angle 
ip). Then, the thermal switching of the magnetization can be 
regarded as the one dimensional Brownian motion of a point 
particle. Although the effect of the demagnetization field of 
an in-plane magnetized system is taken into account in the 
definition of the critical current of Ref. [9], the model of Ref. 
[9] should be regarded as the identical with the models in 
Refs. [10], [11] because the assumption H || p in Ref. [9] is 
valid for the two special cases mentioned above, where H 
and p are the total magnetic field acting on the free layer and 
magnetization direction of the pinned layer, respectively. 

The difficulty to calculate the spin torque assisted thermal 
switching probability arises from the fact that the spin torque 
cannot be expressed as the torque due to the conserved 
energy. Mathematically, it means that we cannot find any 
function F(9,ip) whose two gradients, 8F /dp and OF/ 89, 
simultaneously give the spin torque terms of the Landau- 
Lifshitz-Gilbert equation in (9, ip) coordinate. Then, the steady 
state solution of the Fokker-Planck equation deviates from 
the Boltzmann distribution. However, in the two special cases 
mentioned above, since the magnetization dynamics depends 
on only 9, F can be obtained by integrating the spin torque 
term with respect to 9. Then, the Fokker-Planck equation, 



8W 
~~8t 



ay 



d 



sin 9 89 



+ 



-^appl 

k B T8W 



TJ 

— +H K cos9\sm9W 
a 



MV 89 



(25) 



has a steady state solution of the Boltzmann distribution form, 
W oc exp[-J7(fc B T)]. Here M, V, H appU H K , i7 s (oc I), 
70 = (1 + ct 2 )^', and a are the magnetization, volume of the 
free layer, applied field, uniaxial anisotropy field, strength of 
the spin torque in the unit of the magnetic field, gyromagnetic 
ratio, and the Gilbert damping constant, respectively. F = 
—MH app \V cos 9—(MHkV/2) cos 2 9 is the magnetic energy, 
and & is the effective magnetic energy given by 



-H appi cos 9 



IT j 

— cos 6»- -H K cos 2 9. (26) 
a 2 



The term — (M H s V/a) cos 9 in Eq. d26l corresponds to F 
mentioned above. By using the steady state solution of the 



4 



Fokker-Planck equation, we can calculate the switching prob- 
ability, according to Refs. [8], [10], [11]. 

Koch et al. argued that Brown's formula with the magnetic 
energy F is applicable to the spin torque switching problem 
by replacing a and T with a = a[l + H s /(aH )] and T = 
T/[l+H s /(aH)}, where H = |H| = |if app i+if K cos 0\. These 
replacements arise from the assumption that the directions of 
the spin torque (oc M x (M x p)) and the Landau-Lifshitz 
damping (oc M x (MxH)) are parallel, i.e., H || p. At the 
minimum of the magnetic energy F, T — T/(l — I/I c ), and 
thus, Ref. [9] argued that the exponent of the current term of 
the potential barrier height (oc M H^V / {2k B T)) is unity. How- 
ever, it should be noted that the definition of the the potential 
barrier height requires not only the minimum of the magnetic 
energy F m [ n — F(0) but also its maximum F max = F(8 m ) 
divided by the temperature, where 8 m = cos _1 (— H apv \/ H^). 
We can easily verify that H, and also T, are zero at 9 = 9 m . 
Thus, F mei y i /[k B T(8 m )} is not well defined, and the relation 
argued in Ref. [9] is not satisfied, as shown below: 



F„ 



{F n 



0(i -TO 



k B T 



(27) 



3 T(0 m ) k B T{0) 

The origin of the problem in Ref. [9] is that exp[— i* 1 / (k B T)] 
is not a steady state solution of the Fokker-Planck equation 
(1251 : the steady state solution is exp[— ,^ / {k B T)]. Since the 
effect of the spin torque can be regarded as an additional term 
to the applied field, as shown in Eq. (|26| |. the potential barrier 
height of the spin torque assisted thermal switching is, similar 
to Brown's formula [8], given by 

-^appl 



1 



H s /a 



k B T u V ' H K 

where the thermal stability is defined by Ao 
MH K V/(2k B T). By using the relation 



(28) 



1- 



appl " 



-H s /a 



1- 



H, 



appl 



i- 



a(H K 



and defining the critical current I c by H s /[oi(Hk 
-I/I c , we find that [11] 



5? _ 3? . 
■ y max ■ y min 



1 



H, 



appl 



l 



-#appl)- 

(29) 

- -ffappl)] 



(30) 



k B T u V' H K 

Thus, the exponent of the current term should be two. 

IV. Conclusion 

In conclusion, we studied the spin torque assisted thermal 
switching of the single free layer theoretically. We derived the 
theoretical formulas of the most likely and averaged switching 
currents of the sweep current assisted magnetization reversal, 
and showed that the value of the exponent b in the switching 
rate significantly affects the estimation of the retention time 
of MRAM. We also discussed the difference between the 
theories in Ref. [9] and Refs. [10], [1 1] from the Fokker-Planck 
approach, and showed that the exponent should be two. 
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